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A nonperturbative theory of multiphonon anharmonic decay of strongly ex- 
cited local mode is developed whereby the mode is considered classically and 

^ ■ 
O 

' phonons, quantum mechanically. The decay rate of the mode is expressed 

O ' 

via the negative frequency parts of the phonon pair correlation functions. In 
the case of two-phonon decay the later satisfy the linear integral equations 
while in the case of two- and more-phonon decay they satisfy the nonlinear 
integral equations. As a result, the processes mentioned differently depend on 
the mode amplitude A: two-phonon processes smoothly deminish if A — > oo 
while three- and more-phonon processes are fully switched-off at large ampli- 
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^ ■ tudes and they abruptly switch-on if the amplitude approaches the critical 



value. At that the decay rate gets rather high value (of the order of the mode 
quantum per period). The final stage of the relaxation is well described by 
the perturbation theory. 

I. INTRODUCTION 

The problem of the time evolution of strong vibrational excitations in crystals nowadays 
attracts a remarkable attention. One of the essential results in this field is the demonstration 
of the existence of localized modes (ILM, intrinsic local mode or SLS, self-localized ; 
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in perfect crystals. The results are obtained by applying equations of classical nonlinear 
dynamics. The central point concerning the ILMs is their stability. In the classical limit 
these modes are stable if all their harmonics are out of resonance with the phonon spectrum. 
This condition is fulfilled e.g. for the local modes with high frequency (the basic frequency 
of the mode is above the phonon band). On the contrary, the quantum mechanics allows 
the distinction of the high frequency quanta (due to the processes of creation of two or 
more phonons). Therefore, the question of the stability of the modes mentioned should 
be considered on the basis of the quantum theory. However, the problem mentioned can 
not be solved by applying standard quantum perturbation theory, commonly used for the 
description of quantum transitions since one cannot assume that the perturbation parameter 
is small (note that one of the essential conditions of existence of ILMs in a nonlinear lattice 
is a large amplitude of vibration). 

Recently we proposed a nonperturbative theory of the two-phonon anharmonic decay of 
strong local modes associated with lattice defects^ and an analogous theory of the decay of 
the ILMs in the perfect monatomic chain0. We have found that the decay law is highly 
nonexponential: at definite "critical" amplitudes (typically of the order of 0.2 - 0.4 A) 
the relaxation jumps take place being accompanied by generation of phonon bursts. The 
process considered ini'0 takes into account cubic anharmonicity. The interaction term which 
is responsible for the process is ~ J2ni 2 W^Qtyx^Xiz, where W^J 2 are the parameters of 
the cubic anharmonicity, Q(t) is the time-dependent classical coordinate of the strong local 
mode, x ix and x i2 are the coordinate operators of phonons. Such interaction leads to the 
quadratic (with respect to phonon operators dj, af) time-dependent phonon Hamiltonian. 
The latter can be diagonalized exactly by means of the linear transformation of phonon 
operators 



the terms ~ af result from fast (oscillatory) change in time of the phonon Hamiltoniana. 
Due to these terms the initial zero-point state |0) is not the zeroth state of the transformed 
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(1) 



destruction operator bj, i.e. phonons are created in the lattice. This transformation allows 
one to find energy of generated phonons and, thus the energy loss of the mode for any cubic 
anharmonicity and at any time moment. 

Let us mention that an analogous transformation of the field operators takes place in the 
field of a gravitationally collapsing star (black hole), causing thermal radiation (so-called 
quantum emission of black holes)0-El. Note also an analogy of emission of phonon pairs by 
a strong mode with superconducting transition: in both cases mean values of destruction 
(and creation) operators (of phonons or electrons (holes), respectively) equal to zero but 
mean values of products of two destruction operators (as well as mean values of products of 
two creation operators) differ from zero. 

The described process of phonon generation by a strong mode is allowed if lo\ < 2o>m, 
where ui is the frequency of the local mode, ojm is the top band frequency. If u>i > 2u>m, then 
one needs to account for three- (or more) phonon processes. These processes arise from the 
interaction terms ~ J2k ^iii^*] _ 1 Q k (t)^ii^i2 ■ ■ ■ &in-i with n > 4. The phonon Hamiltonian 
with such terms cannot be diagonalized. Therefore, the method of! is not directly applicable 
for n > 4. Below we present another method which is based on the asymptotical behaviour 
of a pair correlation function and which is applicable for arbitrary n. 



II. GENERAL FORMULAS 

We take into account that according to the Coleman theorem@ the vacuum (zero-point) 
state of bosons in the classical field is not invariant with respect to the transformation 
which changes the classical field. This means that destruction operators are different for 
different values of the classical field. Therefore, if one describes the effect of a classical 
mode, anharmonically interacting with phonons in the quasi-harmonic (Hartree-Fock, or 
mean-field) approximation 

H ph (t) = Y;^{b+{t)b 3 {t) + 1/2) (2) 

j 
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then the destruction operators linearly change with the mode amplitude. For slow relaxation 
(as compared to the mode frequency ui) the transformation is diagonal in the phonon space: 

bi(t) ~ fii{t)ai + i/i(t)af 

(fj,ij ~ HiSij^Vij ~ ViSij] see (1)). Here 6$ and af are phonon operators before excitation 

of the local mode (which is supposed to take place at time t = 0); \fif\ — = 1; the 

constant term on the right hand side of (2) is considered to be small and neglected. In 

this approximation all multiphonon decay processes are described in the same way: the 

relaxation rate of the mode 

dE t {t) = dE ph (t) = ^ d\vj(t)\ 2 
dt dt Y Ui dt 

is determined by the negative frequency terms in the large time asymptotics of the two- 
time correlation function 

< 0\xi(t + r)xi(t)\0 >- ^-[|^(t)| 2 e _iWiT + Wi{t)\ 2 e^ T ], t » r ~ cuf 1 (4) 

(here fast oscillating terms ~ exp (±iuji(2t + r)) are averaged out). 

|z/j(t)| 2 can be found as follows. Presenting coordinate operators by means of initial 



creation and destruction operators £«(£) = yh/2u)i(gi(t)di + g*(t)d^) one gets 

(0\xi(t + r)xi(t)\0) = (h/^Mt + r)g*(t). (5) 

The time dependence of the coordinate operators can be determined from the Lagrangian 
equations of motion. Inserting these equations to the correlation functions (5) and expressing 
multiphonon correlation functions via the Wick-type products of pair correlation functions 
one gets linear (for cubic anharmonicity) or nonlinear (for quartic and higher order anhar- 
monicity) equations for the pair correlation functions. It is straightforward to show that for 
the multiphonon decay one obtains the same equations if to calculate g^t) from classical 
equations of motion for dimensionless normal displacements with " nonclassical" initial con- 
ditions gi(t) = exp {—iujjt + fa), t < and then to average the correlation functions over the 
random phases 4>i (the later averaging is denoted below as (...)). 



We demonstrate the method for anharmonic interactions in the pair potential approxi- 
mation with account of the central forces only. For the case of simplicity we suppose that 
the strong odd mode is fully localized on the site I = 0. Taking Q(t) = A t cos uit one gets 
the following interaction Hamiltonian 
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Htt 1] * ^cosutf n > 2, (6) 



where = = - £ fe W^ k+1) A k t / {2{n - 1)!) and u t = u lt depend on the mode 

amplitude, q m (t) = J2i^miXi(t). Lagrangian equations for g^t) are the following: 

9i (t) + uf gi (t) + 2A t coBUit^UieiraV^g^it) = 

m 

(e m i = e m i/u>l ). Corresponding integral equations read 

9i (t) = e^+tx + 2A t £ e mt V^ f dt x sin(^(t - t^g^h) cosufa (7) 

where g m (t) = J2i^imgi(t) satisfy the analogous equations 

9m(t) = f m (t)+2A t J2 I dhG^it-tjV^g^ih) cos ufa, (8) 

, Jo 



Gmm'it) = J2i ^im^im' sin oj{t is the dynamical Green functional, f m (t) = J2i ei m exp(—iujit + 



Substituting (7) to (5) and taking into account (3) and (4) one gets for t ^> u l 1 



dE ph (t) _ hA 2 t (n - 1] 



dt 2tx Jo 

where 



dujuSpilmGiut - uo)V {n) D^it; u)V { % (9) 



D { :J\t; = [°° dre^D^t; r), (10) 
is the positive frequency (a; > 0) part of the n — 1 - phonon spectrum; 

Anm'(t; r) = (g m (t + r)g* m ,(t)), (11) 
G mra '(w) = / e luJ G mm '(t)dt 



is the phonon Green function in frequency representation. In (9) it is accounted that the main 
contribution to the relaxation rate is given by tx, t[ ~ t and that the negative frequency part 
of the correlator (gi(t+r)g*(t)) comes from the negative frequency term A\ exp(iui(tx — t'x))/^ 
of the Q{tx)Q{ty) ( here t\ and t[ are the integral variables in (7) and (8), t\ > t[). 

III. TWO-PHONON DECAY 

In the case of two-phonon decay (n = 2) the equations (7) are linear. These equations 
have been solved in! which allowed us to derive the explicit expression for the relaxation 
rate. Here we present another derivation of the formula for the relaxation rate which is 
based on the direct calculation of the pair correlation functions D<V(t;w) = D(t;u) and 
which can be extended for three- and more- phonon decay processes. 

Taking into account (8) and (11) one gets the following equation for D(t; r): 

D(t; r) = d*(t; -r) + 2A t f*° dtxG(t - tx)V®D(t; h - t) cos^ti, (12) 

Jo 

where t = t + r,d(t;r) = (g(t + r)f*(t)). For t 3> r < uf the matrix-function d(t;r) 
satisfies the same equation as D(t;r) with d (r) = (/(t)/*(0)) instead of d*(t; — r). This 
equation (as well as (12)) can be solved as follows. First we present this equation in the 
form 

d(t; t) = do(r) + 2A t [° dtxG(t Q - hW^doitx - t) cosu^ + 
Jo 

AA 2 t (° dtx [ 1 dt 2 G{t - ti)y (3) G(ti - t 2 )V {3) d{t; t 2 - t) cosu^ cosufc. (13) 
Jo Jo 

We are interested only in the positive frequency parts of D(t; r) and d(t; r). The linear (with 
respect to A t ) term in the right-hand side of (13) gives only the fast oscillating contribution 
to these parts and it averages out. Omitting also the fast oscillating part of the last two 
terms in (13) one gets the following equation: 

d(t; t) ~ d (r) + A\ r° dtx t dt 2 G(t - tx)V®G(tx - t 2 )V (3) d{t; t 2 - f) e ^^-ti) ( 14 ) 
Jo Jo 

Introducing 
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d m (t;u) = J _ t e^ T d (0) (t;r) (15) 



one obtains 



roc ptg pt\ 

d(t;u) = do(t;u) + A 2 dt dt x dt 2 e lU}{ - t °~ t) 

Jo Jo Jo 

G(to - £i)y (3) G(ti - t 2 )V^d(t; t 2 - t) exp(i^(t 2 - *i)) (16) 
Denoting to = to — t\, ii = t\ — t 2 and taking into account that 

roc rto rt\ roc roc roc 

/ dt / dt\ I dt 2 = dt 2 / dti I di 
Jo Jo Jo Jo Jo Jo 

one gets simple algebraic equation for d(t; u) with the following solution: 

d(t; w) = [I- J^G(-u)V^G(ui - u)V®]- ld ofc (17) 
Equation (12) can be solved in the same way. As a result one finds 

D (1 \t;uj) = -\I - AlGi-^V^G^- uo)V i3) \~ 2 ImG(uj) (18) 
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and 

™*® = .*m . r d^iv^imG^ - U ) 

dt dt 2tt Jo y[ K ' 

V {3) \I - A 2 G(uj - uoi)V {3) G{uo)V {3) \- 2 ImG{uo)}. (19) 

Within notations and cyclic permutation of the matrices under Sp this formula coincides 
with analogous formula derived in! by another method. Note that diagonal elements of the 
resolvent-matrix [I — A 2 G(uo — io{)V^ G{lo)V^ ) \~ 1 are symmetric with respect to the change 
uj <r-> uj — Taking also into account that G(— u) = G*(u) one finds that the resolvent 
([I - A 2 G(-ui/2)V^G(ui/2)V^]- 1 ) mm is real. Therefore, the function under integral (19) 
for some critical values of A t diverges at u = ui/2. This leads to sharp peaks of the relaxation 
rate near A cr of the shape \A t — A cr | _1 ~ \t — t cr \~ l l 2 where time moment t cr of the jump 
corresponds to A t = A cr i. The relaxation jumps are accompanied by explosion-like emission 
of phonon bursts. 



IV. THREE- AND MORE-PHONON DECAY 



Let us consider now decay processes with creation of more than two phonons, caused 
by quartic and higher order anharmonicities. Substituting (8) to (11) one obtains after 
averaging over ^ the following nonlinear integral equations for the matrix-function D(t; r): 

D(t; r) = d (r) + 4(n - 1) f° dt x f dt 2 Q(t 1 )Q(t 2 )G(t - tjV™ D^fa h - t 2 )V^G(t - t 2 ), 

Jo Jo 

(20) 

where Z)( n_1 ) (t; t) is matrix with the elements (D^ n ^(t; r)) mm / = D r ^j{t]T) (in (20) it is 
accounted that terms ~ {f(t + T)(g n ' 1 (t'))*) and (g^ 1 ^')/* (*)) turn to zero for n > 3 when 
averaging over random phases (pi). The Green functions G(t) essentially differ from zero 
only for small |r| < uf 1 . Therefore, in the limit t ^> r < ujf 1 one can replace the lower 
limits of integration by — oo. Neglecting also fast oscillating terms ~ cos(c<^(ti +t 2 )) one 
gets 

D(t;r) ~ do(r) + 2(n- 1)A 2 / / cfocos (wj(r - a;)) 

JO J-Tl 

r — x)l / ^ n ' ) G'(r 1 ) (21) 
(here Ti — t — t 2 , x — t — t\ — t±). Taking into account that 

roo roo roo roo rO rco 

/ dri I dx — l dx l dr\ + dx dri (22) 

Jo J-ti Jo Jo J-oo J-x 

and introducing the fourth-rank matrix-function 

(x) — / dTiG mmi {ji + x)G m i m2 {ji) (23) 
Jo 

one gets the following set of nonlinear integral equations: 

D mm >(t', T) = G? , mm '(T) + 

2(n-l)A 2 t £ V<pV$> f~ dxca&(unx)D^frx)G^^(\T-x\) (24) 

(n > 3). These equations can be solved numerically and, in some cases, analytically. Then, 
performing Fourier transformation (10) one can find D( n ~ 1 )(t; to) and then the relaxation 
rate (9). 



V. SWITCH-OFF OF MULTIPHONON PROCESSES 



As an example, let us consider a n-phonon decay (n > 3) of a high frequency local vibra- 
tion, directed to the nearest neighbours. In this case only two coordinates q m of the nearest 
neighbours can be accounted. Besides, the high frequency local mode breaks the dynamical 
correlation (nonorthogonality) of these coordinates. Due to that the Green functions G(t) 
and become diagonal. Therefore the equation (24) takes the form 

/OO 
dxD n ~\x)G(\r-x\) cos oj iX (25) 
-OO 

(index m is omitted). We consider the case of phonon spectrum with narrow single maximum 
at ujq = ijJijn of the Lorentzian shape with the width T(^C cu ). 
In this case 

G(t) = wo 1 e _r|T| sin u r; d (r) = c^V^-^ 1 ; G(x) = -^e~ r ^ cos (uqx) (26) 
Approximate solutions of the integral equation (25) is 

D{t; t) « uo^ate-^ + Ae^ or ]e" r|r| (27) 

(t > t ~ ujq 1 ), where f3 t = w t a r ^ 1 , 

h n - 2 V^ 2 A 2 t (n-l) 2 

v>t= nr2 ( tK , ' ; n>3 28 
2uq T 2 n(n — 2) 

a t is a solution of the equation 

a t = l + <a t 2(n_1) . (29) 



TABLE 1 



n 


3 


4 


5 


6 


7 


8 


9 


10 




0.32 


0.405 


0.47 


0.52 


0.56 


0.595 


0.62 


0.64 


®-n,cr 


1.25 


1.173 


1.127 


1.09 


1.068 


1.059 


1.042 


1.031 
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This equation has real solutions only for w t being below the critical value u> njCr and, corre- 
spondingly for A t < A ncr and a t < a„ iCr ; values of w n)Cr and Table 1). 

There are two such solutions but only the smaller one corresponds to slow relaxation (as 
compared to u>i) and, therefore only this solution should be accounted. The relaxation rate, 
then equals 

(h^r-ir =hUlTa t Wt l—i\- (30) 

For small amplitude A t , 7 = d\nE p h(t)dt does not depend on A t , i.e. relaxation is expo- 
nential in accordance with perturbation theory0. For A t < A n>cr the relaxation rate grows 
as A t — > A nfir and gets rather high value ~ TTioji (of the order of the mode quantum per 
period). Abscence of real solutions of the equation (29) for w t > w nfir (A t > A^\ T ) means 
that n-phonon relaxation is fully switched-off. Note that in the case of two-phonon decay the 
switch-off effect also exists but it takes place only asymptotically with A t — > ooi. The values 
of A n , cr ~ h^^w^lTu^/lV^l for n > 3 are much larger than the corresponding values 
of the critical amplitudes for two-phonon decay (for n > 3 A n>cr — > 00 if h — *■ 0) . Typically 
^3,cr ~ lA -C A^ cr v45 jCr etc. Note also that the the difference of the perturbative and 
nonperturbative treatments on the n-phonon decay for n > 3 is given by the factor a" -1 
in (30). For all A t < A n ^ cr this factor does not exceed 1.62. Therefore, the main results 
of the nonperturbative consideration of the multiphonon decay are the switching-off of the 
process if the mode amplitude exceeds the corresponding critical value, and the step-vise 
switching-on of the very fast relaxation when the amplitude reaches this value from above. 
For not very large amplitudes (< 1A) three- and more- phonon decay processes are well 
described by the standard perturbation theory. 

Thus, the general scenario of the relaxation of strongly excited local mode is the following. 
In the beginning the relaxation (which takes place due to creation of n > 4 phonons) is rather 
slow. Then, when the mode amplitude approaches the critical value for the (n — l)-phonon 
decay processes, relaxation is jump-like enhanced (supposing these processes are allowed 
by the energy conservation law). After that the (n — 1)- phonon processes will slow-down 
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till the next critical amplitude (which coresponds to the (n — 2)-phonon processes), then 
the relaxation process will be once more jump-like enhanced (supposing the (n — 2)-phonon 
process is allowed) and so on. The relaxation rate can change very fast not only for a small 
change of the mode amplitude but also for a small change of the mode frequency. Final 
stage of the anharmonic relaxation is exponential and is well described by the perturbation 
theory. 

In conclusion, we developed here a nonperturbative theory of multiphonon anharmonic 
decay of a strong local mode. We found that three- and more- phonon decay processes of the 
mode are described by essentially different equations than the two-phonon decay process. 
The main effect which is predicted by the nonperturbative theory is the switching-off of the 
n > 3— phonon process if the mode amplitude exceeds the corresponding critical value and 
sharp (step-wise) switching-on of the fast relaxation when the amplitude reaches this value 
from above. 

The author is grateful to A.Bussman-Holder for invitation to participate in the confer- 
ence. The research was supported by Estonian SF Grant No. 2274. 
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